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We determine the contributions from the direct Coulomb and exchange interactions to the total 
interaction in semiconductor artificial atoms. We tune the relative strengths of the two interactions 
and measure them as a function of the number of confined electrons. We find that electrons tend to 
have parallel spins when they occupy nearly degenerate single-particle states. We use a magnetic field 
to adjust the single-particle state degeneracy, and find that the spin-configurations in an arbitrary 
magnetic field are well explained in terms of two-electron singlet and triplet states. 
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The addition of a single electron charge to a quan- 
tum box costs a certain energy, which is responsible 
for Coulomb blockade in electron transport [1]. Also a 
change in spin is associated with a certain change in en- 
ergy, e.g. exchange energy is gained when electrons are 
added with parallel spins as compared to anti-parallel 
spins. Depending on the system, a large total spin (fer- 
romagnetic filling) or a minimum total spin value (anti- 
ferromagnetic filling) is favored. In semiconductor quan- 
tum dots alternate spin filling [2] as well as spin-polarized 
filling [3] have been reported. Here, we study verti- 
cal quantum dots which have well-defined single-particle 
states. When these states are separated by a large en- 
ergy, A_E, an anti-ferromagnetic filling is favored. For 
small /S.E, a ferromagnetic filling is observed, which is 
in line with Hund's first rule from atomic physics. We 
use a magnetic field, B, to tune AE{B) allowing us to 
alter the spin filling. We first discuss a simple model 
that describes filling of two single-particle states with 
two interacting electrons. Fig. 1(b) shows two, spin- 
degenerate single-particle states with energies Ea and Eh 
crossing each other a.t B — Bq. The ground state (GS) 
energy, U{1), for one electron occupying these states, 
equals Ea for B < Bq and Ei, for B > Bq (thick line 
in Fig. 1(b)). For two electrons we can distinguish 
four possible configurations with either total spin S — 
(spin-singlet) oi S = 1 (spin-triplet). (We neglect the 
Zeeman energy difference between Sz ~ ~ 1 , and 1 . ) 
The corresponding energies, Ui{2, S) for i = 1 to 4, are 
given by: Ui{2,0) = 2Ea + C„a,C/2(2,0) = 2Eb -f Ctb, 

[73(2, 1) = Ea+Eb-^Cab-lKabl Ui{2,Q) = Ea +Eb+Cab. 

Here, Cij = a,b) is the direct Coulomb (DC) en- 
ergy between two electrons occupying states with ener- 
gies Ei and Ej, and Kab is the exchange (EX) energy 
{Kab < 0) between two electrons occupying Ea and Eb 
with parallel spins [4]. The experiments below measure 
the electrochemical potential defined for a two electron 
system as ^(2) = U{2) - U{1). For each C/,(2) we ob- 
tain the potentials: Hi{2) — Ui{2) — Ea for B < Bq and 



fi^{2) = U.,{2) - Ebior B > Bo (see Fig. 1(c)). 
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FIG. 1. (a) Scanning electron micrograph of the semicon- 
ductor quantum dot device. The quantum dot is located 
inside the 0.5 /^m-diameter pillar. The top and bottom 
contacts serve as source and drain electrodes. A Schottky 
gate is wrapped around the pillar. The current, 7, fiow- 
ing through the dot is measured as a function of gate volt- 
age, Vg, in response to a dc voltage, V, applied between 
the source and drain, (b) Schematic diagram of two sin- 
gle-particle states with energies Ea and Ei, crossing each other 
at a magnetic field B — Bq. (c) Electrochemical potential, 
/ii(2) = Ui(2) — U{1), for two interacting electrons. The thick 
line depicts the ground state energy whereas the thin lines 
show the excited states. 

The GS has S — away from Bq. Near Bq, the lowest 
energy is /is (2), such that here S* = 1. The downward 
cusp in the thick line identifies this spin-triplet region. 
The transition in the GS from S* = to 1 and 5 = 1 to 
0, respectively, occurs when fil = for B < Bq and 
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when /i2 = /Z3 for B > Bq. We define two energies, Ai 
and A2, to characterize the size of the downward cusp in 

the GS at B = Bq: Ai = - /ig = Caa ~ Cat + \Kab\, 
A2 = — A*3 = Cbb — Cab + \Kab\, Ai — A2 = Caa — Cbb- 

Our semiconductor quantum dot (see Fig. 1(a)) has the 
shape of a two-dimensional disk [5]. For detecting the 
GSs, we set the source-drain vohage, V, to a small value. 
The GSs and excited states (ESs) are both measured 
when V is set to a value sufficiently greater than the ex- 
citation energy [6]. For small V ^ a series of current peaks 
results from changing the number of electrons in the dot, 
N, one-by-one [1]. The position of a current peak for 
the transition from iV — 1 to iV measures the GS electro- 
chemical potential fJ.{N). The sample is cooled down to 
about 100 mK. Fig. 2(a) shows the evolution of current 
peaks with magnetic field for = 7 to 16. 
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FIG. 2. (a) Evolution of the ground state energies from 
A'' = 7 to 16 as measured from the current peaks versus mag- 
netic field at V = 120 mV. The bars along the gate voltage 
axis show 1 meV energy scales calibrated at -1.26 and -0.85 
V. The dotted curve indicates the last crossing between sin- 
gle-particle states. Dashed ovals correlate pairs of ground 
states for odd and even electron numbers. Spin transitions in 
the ground states are indicated by ■ at B = T and occur in 
the ovals for B 7^ T. (b) Magnified plots of the A'' = 8 and 
24 current peaks vs. magnetic field. The dashed lines illus- 
trate how the interaction-energy parameters, Ai and A2 are 
determined, (c) Fock-Darwin single-particle states calculated 
for ujo — 2 meV. 

Features associated with a parabolic confining poten- 
tial are all observed such as a shell structure [2]. The 
large spacing for iV = 12 at _B = T can be seen in Fig. 
2(a) (see double arrow) and marks the complete filling 
of the first three shells. The pairing between neighbor- 
ing peaks indicates anti-parallel spin filling of a single 
orbital state by two electrons. Modifications to this pair- 



ing are observed for the peaks labeled by ■ at T, and 

in each of the dashed ovals connecting pairs of peaks at 
non-zero field. These are all signatures of Hund's first 
rule; i.e. spin-polarized filling. Note that the Zeeman 
effect is negligible in this experiment [7]. We show ex- 
pansions of the evolution of the = 8 and A^ = 24 
peaks in Fig. 2(b). The downward cusps are clearly 
seen. The dashed lines form a parallelogram, from which 
we obtain parameters Ai and A2. To compare the two- 
electron model with larger electron numbers, we assume 
that other states are far away in energy so that they can 
be neglected. Then, the downward cusps should occur 
for higher euen-electron numbers, whereas they should 
be absent for odd-electron numbers. This is clearly ob- 
served in the ovals in Fig. 2(a). For instance, the B-field 
dependence of the 9th peak compares well to the thick 
line in Fig. 1(b) and the B- field dependence of the 10th 
peak compares well to the thick line in Fig. 1(c). Other 
pairs of even and odd numbered peaks show the same 
behavior. 




FIG. 3. dl/dVg in the plane of Vg and B for A^ = 7 
to 9 measured for V — 2 mV. dl/dVg > for dark blue, 
dl/dVg < for white and dl/dVg ~ for yellow. The solid 
red lines indicate the evolution of the GSs with magnetic field 
whereas the dashed red lines show the ESs. The two arrows 
indicate singlet-triplet (S-T) and triplet-singlet (T-S) transi- 
tions in the GS for A' = 8. 



This justifies our assumption so that we can simplify 
the many-electron system to just one or two electrons. 
More detailed agreement is obtained by measuring the 
excitation spectrum [6]. Fig. 3 shows a dl/dVg plot, 
taken for V — 2 mV. This larger voltage opens a suffi- 
ciently wide transport window between the Fermi levels 
of the source and drain, that both the GS and first few 
ESs can be detected. The GS and ESs for TV = 7 to 9 
can be assigned from the magnetic field dependence of 
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the dark blue lines. Solid red lines highlight the GSs 
whereas the ESs are indicated by dashed red lines. TThe 
set of GS and ES lines for TV = 7 shows a single crossing 
similar to that in Fig. 1(b). The spectrum for N = 8 
compares well to Fig. 1(c) and we can clearly distinguish 
the parallelogram formed by the GS and first ES. The 
downward cusp in the GS for N = 8 (labeled A) is at a 
slightly higher B-field than the upward cusp in the first 
ES (labeled ▼). This asymmetry implies that Ai > A2, 
i.e. Caa > Cbb. The same type of asymmetry is always 
observed along the dashed line in Fig. 2(a), implying that 
Caa > Cbb for all N. Note that the GS for TV = 9 shows 
an upward cusp (labeled by V) quite similar in form to 
the first ES in the spectrum for = 8. This implies 
that the filling of the ninth electron is closely linked to 
the configuration of the N = 8 first ES [3] . As illustrated 
in Fig. 2(b), we can derive the experimental values for 
Ai and A2 for different N. These values are plotted in 
Fig. 4. We find that Ai is larger than A2 for all N, again 
implying that Caa > Cbb. As N increases from 6 to 12, 
Ai first increases and then slowly decreases, whilst A2 
slightly decreases. 
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FIG. 4. Experimental values for the energy parameters Ai 
(•) and A2 (■), on a log. scale, versus electron number de- 
rived from data as shown in Fig. 2(a). The uncertainty in 
the determination of the eperimental values is lOcalculated 
from the FD-wavefunctions at 2 T, for an unscreened (unscr) 
and screened (scr) Coulomb interaction The calculated ex- 
change energy li^atl between states with energies Ea = -Bo,-i 
and El, — -Eo_jv/2-i decreases quickly with A'^. Above the main 
figure the absolute squares of the wavefunctions are shown for 
the relevant quantum numbers n — and / — 0, ±1, ±2, ±3. 
As the angular momentum quantum number I increases, the 
average radius increases. 

To calculate DC and EX energies [4] we now need to 
specify the confining potential of the disk. In earlier 



work it was shown that the lateral confinement is well 
described by a parabolic potential with cylindrical sym- 
metry [2]. The eigenfunctions with eigenenergies, 
in this potential are known as Fock-Darwin (FD) states 
[8]: 

E,u^-^-hu, + {n + ^ + ^\l\)h,/i^^!l+Zj (1) 

where rt = 0, 1,2, ... is the radial quantum number and 
I = 0, ±1, ±2, ... is the quantum number for angular mo- 
mentum, hiuo is the lateral confining energy and HlOc 
= eB/m* is the cyclotron energy. Each FD-state is spin 
degenerate. At i? = T the FD-spectrum has sets of 
states with increasing degeneracy (see Fig. 2(c)). This 
degeneracy is lifted on increasing B, but as B is increased 
further new crossings can occur. The last crossing is al- 
ways a crossing between just two FD-states. The up- 
going state is always (n, I) — (0, —1), whereas the down- 
going state, (0, 1 > 1), has an increasing angular momen- 
tum for states with increasing energy. (The relation with 
Fig. 1 is: Ea = Eo,-i and Eb = Eo,i>i ) Note that the 
last crossings also correspond to the dashed line in Fig. 
2(a). From the electron distributions of the FD-states 
we calculate the DC and EX energies for two electrons 
occupying two degenerate states. We take hujo = 2 meV 
as deduced from earlier experiments [2,6] and obtain Ai 
and A2. The dashed curves in Fig. 4 show Ai and A2 
when we neglect screening of the interactions within the 
dot by electrons in the leads and in the gate. In this case 
the Coulomb potential falls off as 1/r, where r is the 
distance between the electrons [4]. For the solid curves 
we have approximated the screening effects by replacing 
the Coulomb potential by exp{— r/dj/r. We have taken 
d = 10 nm which is roughly the thickness of the tunnel 
barriers. Fig. 4 shows that screening considerably re- 
duces Ai to values much closer to the experimental val- 
ues. Screening also removes the minimum in A2, which is 
also in better agreement with the experiment. Since the 
average radius of the wavefunctions increases with angu- 
lar momentum, two electrons are closer together when 
they both occupy (0,-1) compared to when they both oc- 
cupy {0,1 = A^/2 - 1) for even N > 4 {or I > 1), so 
the DC interaction is stronger in the former. This ex- 
plains our observation (Ai > A2) Caa > Cbb for all N. 
The overlap between different wavefunctions, (0,-1) and 
(0, ; = N/2 - 1), decreases for even N > i {or I > 1). 
This results in a decrease in both Cab and \Kab\ with N. 
It then follows that Ai increases until it saturates at a 
value equal to Caa- The gradual decrease of experimen- 
tal Ai for iV > 12 is probably related to the decrease 
in the lateral confinement with N [2] and thus the de- 
crease in Caa- We finally discuss the interaction effects 
for the iV = 4, 8 and 14 peaks near B = 0T (the iV = 8 
and 14 peaks are labeled n in Fig. 2(a)). These cor- 
respond to the GS electrochemical potentials for adding 
the second electron to the second, third and fourth shells. 



3 



respectively. The inset to Fig. 5 demonstrates the resem- 
blance to the model of Fig. 1(c) for N = 8 near B = 
T. Comparing these data to the FD-spectrum, we assign 
the states such that: Ea = -E0.-2 and Ei, = i?o.2- Like- 
wise, for = 4 we have Ea = i?o.-i and Ei, = £^0,1 [2] 
and for = 14 we have Ea = i?o,-3 and Eb = £'0,3 [9]. 
Note that these states correspond to wavefunctions with 
a complete overlap. Also, for B = Bq = T the two 
crossing states have the same orbital symmetry implying 
Ai = A2 = |ifoh|; i-e. only EX effects contribute to the 
downward cusp [10]. 
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FIG. 5. Exchange energy [ifabl, on a log. scale, associ- 
ated with spin-triplets formed when each new shell is filled by 
just two electrons at B = T. The solid circles are the ex- 
perimental values whose uncertainty in the determination is 
lOthe filling of the first two electrons into the third shell (i.e. 
N = 7 and 8). The vertical double-arrow represents \Kab\ in 
units of gate voltage which is then converted to energy. The 
calculated curves in the main figure are for the unscreened 
(dashed) and screened (solid) cases. 

We derive ji^ablas illustrated in the inset to Fig. 5. 
The obtained EX energy quickly becomes smaller for 
higher lying shells. For comparison we also show the 
calculated screened and unscreened values. The screened 
case provides the best quantitative agreement for our re- 
alistic choices of the confining energy and the screening 
distance. Our general model provides a clear identifi- 
cation of effects due to EX and DC interactions. More 
advanced calculations support our analyses [11]. An im- 
portant simplification is the reduction of a many-electron 
system to just two interacting electrons. The type of spin 
filling in many nearly degenerate levels near i? = T in 
larger electron boxes remains an interesting open issue. 
We thank G. Bauer, M. Danoesastro, M. Eto, R. van der 
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